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Recent biochetnical and clinical evidence implicates human serum camosinase in a variety of pathological 
conditions, such as neurological disorders and diabetic nephropathy, suggesting that this enzyme is of potential 
mterest as a novel medicinal target. The present study was undertaken with a view to model the serum 
camosinase and its catalytic site and to unravel the molecular mechanism by which citrate ions increase the 
catalytic efficiency of seram camosinase. A homology model of the enzyme was obtained on the basis of 
^-alanine synthetase, and its active center was found to bind known substrates camosine, homocamosine, 
and ajisenne in a binding mode conducive to catalysis. Citrate ions were shown to bind at only three weU- 
defined sites involving both ion paks and hydrogen bonds. MoJecular dynamics simulations evidenced that 
citrate binding had a remarkable conformational influence on the 3D structure of camosinase, increasing 
the binding affinity (i.e., binding score) of camosine to the catalytic site. This is one of the first reports 
documenting the molecular mechanism of an allosteric enzyme activator using MD simulations. 

Introduction 

The dipej^tide camosine (^-alanine-L-histidine) represents the 
archetype of a series of histidine-containing dipeptides, such 
* as homocamosine, carcinine, TV^-acetylcamosine, and anserine 
(Chart 1).^'^'^ The complete role of these dipeptides is still 
unknown, even though their function has been studied inten- 
sively in recent years. Available studies indicate that camosine 
can act as a cytosolic buffer,"^ an antioxidant^^ an a,^-.unsatured 
carbonyl scavenger,^ and an antiglycation agent"^ Moreover, 
homocamosine can represent a GAB A reservoir^ mediating the 
antiseizure effects of GABAergic therapies.^ 

Camosine is synthesized by camosine synthetase^ (EC 
6.3.2.11) in many tissues^^ and hydrolyzed by ubiquitous 
dipeptidases (also called camosinases), which belong to met- 
alloproteases.^^ Until now, two isofomis of camosmase have 
been characterized. The first enzyme is a cytosolic form (also 
named tissue camosinase, CN2, EC 3.4.13.18), which acts as a 
nonspecific dipeptidase with broad substrate specificity and 
being strongly inhibited by bestatin.^^ jj^^ second enzyme is 
known as serum camosinase^^ (CNl, EC 3.4.13.20) and is 
characterized by its distribution in the plasma and brain, its 
abihty to also hydrolyze anserine and homocamosine, and its 
absence in nonprimate mammals except in the Syrian golden 
hamster.'"* Recently, the cDNA codes for the two carnosinases 
were identified and expressed, confirming many biophysical and 
biochemical properties of the two isoforms.*^ 

Experimental evidence shows that senim camosinase is 
present in solution as a homodimer, despite the fact that the 
monomer alone appears to be catalytically self-sufficient 
Although tlie nature of the metal ion in serum camosinase 
remains unknown, it is likely tliat two Zn-*^ ions are involved 
in the catalytic site because of its homology with die M20 family 
of metalloproteases.^^ Finally, semm camosinase is activated 
by Cd^"** and citrate ions,'^ although it is unknown whether 
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Chart 1, Histidine-Containing Dipeptides Investigated Here 
camosine \ss^N 




camosine homocamosine 



^^^N carcinine ^^v^^'^!^ 



N-acety!camoslne 



^CHg * • anserine 

citrate promotes enzyme dimerizadon and/or improves the 
catalytic ej65ciency of the single monomer by acting as an 
allosteric activator. 

There is a large body of evidence demonstrating that serum 
camosinase is involved in some important pathological condi- 
tions. Thus, low levels of serum camosinase is responsible for 
elevated concentrations of homocamosine in the brain and 
cerebrospinal fluid (homocamosinosis),^^ a deficiency ac- 
companied by neurological disorders, such as paraplegia, retinitis 
pigmentosa, and progressive mental deficiency.^'^ Decreased 
concendrations of semm camosinase have been also observed 
in patients with Parkinson's disease, multiple sclerosis, or after 
a cerebrovascular accident.'^ It has also been suggested that the 
monitoring of serum camosinase may be useful to predict the 
clinical outcome of patients with acute stroke. Furthermore, 
a recent study^*^ has demonstrated that diabetic patients with 
genetically based low serum camosinase levels {Mannheim 
variant) are less susceptible to nephropathy,- suggesting some 
beneficial effects of carnosine against the adverse effects of high 
glucose levels in renal cells. 
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Figure L Three-dimensional ribbon of the two models of human serum camosinase generated in this study. The lid domain is colored in green, 
whereas the catalytic domain is in blue. The metal ions are displayed in silver. (A) Three-duDensional model based on jS^-alanine synthetase. (B) 
Three-dimensional model based on PepV. 



Taken globally, these studies shed light on the proraismg roles 
of serum camosinase and suggest a therapeutic useftilness of 
either inhibiting (e.g., in diabetes) or activating (e.g., in homo- 
camosinosis) the enzyme. The design of seiTim carnosinase in- 
hibitors can be based on our detailed knowledge of its catalytic 
site and its recognition capacities. In contrast, almost no infor- 
mation is currently available to design semm camosinase acti- 
vators, but a possible path to such information may be opened 
by exploring the molecular mechanism by wliicb citrate activates 
the enzyme. Indeed, such an understanding may unveil novel 
approaches to increase the activity of serum carnosinase, for 
ejcample, by permitting the rational design of citrate-like, non- 
toxic allosteric modulators. 

To approach this objective, we first undertook to generate a 
3D model of human serum camosinase by homology modeling 
techniques. The resulting model was then validated by docking 
a few histidine-containing dipeptides, thus unravelling the 
phamiacophoric pattern of the catalytic site. In a third step, 
molecular dynamics (MD) simulations were used to exanme, 
ill molecular detail, the effects of citi-ate ions on the activity of 
semm camosinase. 



Results 

Structural Features of Camosinase. As described in the 
Computational Methods section, two homology models of serum 
camosinase were generated using PepV dipeptidase and y3-ala- 
nine synthetase as templates. However, the first model was 
discarded because of the inaccessibility of its catalytic site, and 
all'docMng analyses and MD simulations were carried out using 
the model derived from ^S-alanine synthetase. 

Figure 1 A depicts die serum camosinase model obtained with 
^-alanine synthetase as the template, showing an ellipsoidal 
shape with axes equal to 80 and 40 A. The 3D structure of the 
protein consists of two distinct domains, namely, the upper or 
lid domain and the lower (catalytic) domain. The latter con- 
sists of residues from Serl to Pro275 and from His390 to 
His481, whereas the lid domain comprises die residues from 
Leu276 to Leu389. Zinc ions are located in the catalytic do- 
main and are involved in interactions with Hisl06, Glyll5, 
Aspl39,- Glul73, Glul74, and two water molecules, as seen in 
/3-alanine synthetase. The water molecules were kept during the 
modeling of the protein but were deleted in docking calculations 
and MD simulations to facilitate substrate interactions. With- 
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of first 10 


ofiirst20 


of all 30 


dipeptide 


hydrolysis*' 


score^ 


frames^ 


frames^ 


frames* 


camosine 


HYD 


-40.25 


37.21 


3539 


34.12 


homocaniosine 


HYD 


"40,21 


37.19 


35.74 


34.41 


carcmine 


■NON-hyd 


-35.08 


32.13 


30.37 


29.12 


iV-acetylcamosine NON-hyd 


-34.10 


29.26 


2735 


25.66 


anserine 


HYD 


""30.75 


25.96 


23.60 


21.99 



" Because of large discrepancies among experimcnial literature data, a 
simple dichotomic classification is used here. HYD = hydrolyzed dipeptides; 
NON-hyd = unhydrolyzed dipeptides. * Score values m kcal/mol. 
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Figure 2. TW(HlimeiisionaI representation of the mteraction pattern 
between camosinase and its substrate camosine. The model shows how 
the enzyme recognizes (binds) tlie ammonium group, the carboxylate 
^up» and the unsubstituted imidazole ting. The anaido bond is 
simultaneously bound for recognition and polarized for catalysis. 

out water molecules, tfie side chain of Aspll6 bridges the two 
zinc ions. 

The catalytic domain consists of a series of antiparallel 
iS-sheets surrounded by a crown of a-helices, whereas the lid 
domain shows some antiparallel /5-sheets in flie interface with 

the catalytic domain and two a-helices in the upper surface. 
Lid and catalytic domains are connected by a random coil region 
that can act as a hinge pennitting the relative movements of 
the domains during substrate recognition and enzyme activation. 
Moreover, the contacting surfaces of the two domains are 
noticeably rich in ionized residues; this generates a number of 
key ion pairs (e.g., Glul87 with Lys346 and Aspl89 with 
Arg350) that stabilize the proper orientation of the domains. 

Figure IB shows that the second serum camosinase model 
based on PepV dipeptidase^^ is comparable to the previous one 
but widi the two domains clearly closer. In particular, this model 
shows a lid domain markedly larger than that in the j3-alanine 
synthetase model consisting of antiparaOel /9-sheets in the 
interface with the catalytic domain and four a-heiices in the 
upper surface. 

Ligand Docking and Phannacophoric Pattern. Table 1 
reports the scores of relevant histidine-containing dipeptides 
when docked with the serum camosinase model based on 
P-^mint synthetase as the template. As seen in Table 1, 
camosine shows the best score among these dipeptides. Interest- 
ingly, the same rank order is obtained considering the best 
solution, the average score of the first 10 complexes, the average 
score of the first 20 complexes, and the average score of all 30 
complexes. 

The best docking mode of camosine was then used to define 
a recognition pharmacophore of serum camosinase toward its 
substrates (Figure 2). The ammonium group of camosine is seen 
to interact with the carboxylates of Aspll6 and Glu451 and 
with the backbone carbonyl of GlyllS. The carbqnyl and 
carboxylate groups of camosine also interact with the Zn^"^ ions 
in the catalytic site. Foiarizadon of the carbonyl group is thus 
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increased, an essential catalytic feature of hydrolases.^ The 
backbone carbonyl of Leu254 forms a hydrogen bond with 
in the imidazole ring. Like Thr424, it produces an interesting 
set of interactions, namely, a H bond with N^, a second H bond 
with the ami do group, and a hydrophobic interaction with the 
imidazole ring. 

It is of relevance to point out that the interactions unraveled 
here involve both nitrogen atoms tn the imidazole moiety, a 
fact that can explain the specificity of serum camosinase toward 
histidine-containing dipeptides. Moreover, these interactions can 
explain the lower activity of serum camosinase toward camosine * 
analogues widi a substituted nitrogen in the hnidazoJe ring, 
aldioogh the capacity of anserine to be hydrolyzcd suggests that 
all imidazole-involving contacts are not required for enzyme 
recognition and catalysis. 

Table 1 shows that bomocamosine has a high dockiug score 
identical with that of camosine, and its mteraction pattern is 
also very similar (data not shown). This result can be justified 
considering that the larger steric hindrance of bomocamosine 
can be compensated by its greater flexibOity. As a result, 
bomocamosine can assume a bioactive conformation similar to 
that of camosine in its best complex, and, thus, be a good 
substrate of the enzyme. This result is of significance because 
it confinns the capacity of serum camosinase to hydrolyze 
bomocamosine and, thus, play an exploitable role In homocar- 
nosinosis-induced disorders (see below). 

Finally, the modest docldng scores of carcinine and N- 
acelylcamosine can be rationalized considering that carcinine 
lacks a carboxylate group, whereas TV-acetylcamosine lacks a 
charged ammonium group and is charactearized by marked steric 
hindrance, which prevents its carbonyl group from interacting 
with the zinc ions. 

Taken globally, die docking results afford an encouraging 
validation of our semm camosinase model based on ^-alanine 
synthetase, and tliey suggest a comprehensive and coherent 
interaction pattern i3etween the enzyme and its ligands. Nev- 
ertheless, it must be remembered diat docking calculations can 
simulate only the recognition phase between enzymes and their 
substrates, whereas the catalytic phase cannot be simulated by 
this approach. In other words, a given compound can be found 
to elicit good binding interactions and high score values with 
an enzyme, yet fail to possess an appropriate target group or 
fail to allow its binding in a catalyticaliy productive mode, as 
seen with //-acetylcamosine (results not shown). 

It is worth emphasizing that the docking calculations were 
also used to select the best carnosinase model. Indeed, docking 
calculations were also performed for camosine using the model 
based on PepV dipeptidase. The complexes obtained with the 
PepV dipeptidase model were of lower quality than those 
obtained with the ^-alanine syndietase model, mainly because 
of the scarce accessibility of the catalytic site m the former 
model (results not shown). 

Docking of Citrate Ions. Docking calculations of citrate ions 
to serum carnosinase led to some unexpected results (Figure 
3). Indeed, the binding of the citrate ion was found to be 
restricted to three high-affinity sites on the protein surface 
(labeled a, fi, and 7 for clarity); all 30 complexes generated 
during the simulations exclusively Involved these diree sites. 
We found it highly surprising tliat the recognition of a small, 
flexible, and strongly negatively charged molecule such as the 
citrate ion by a protein surface rich in positively charged residues 
should occur exclusively on a few well-defined sites. 

Figure 3 shows the best docking score for each binding site, 
revealing that tlic a site is located in the region connecting the 
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Figure 3. Solution clusters obtained with the 30 complexes between 
camosinase and the citrate ion. The Figure shows the best solution for 
each binding site. 
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Figure 4. RMSD profiles for the four MO simulations perfonned with 
catnosinase (CNl), camosinase with camosine, camosinase with three 
citrate ions, and camosinase with camosine plus three citrate ions. 

two domains of serum camosinase, whereas the ^ and y sites 
are located in the lower region of the catalytic domain. In all 
cases, stabilization of the complexes was through io0 pairs with 
positively charged residues (e.g., Lys405 and Arg409 in the a 
site, Lys28 and Arg50 in the p site, and Arg65 in the y) and 
reinforced by H bonds (e.g., Tyr401 in the a site, Gly25 and 
Tyrl23 m the site, and Gly63 and Ser94 in the y site). Indeed, 
the citrate ion never elicited three ionic bonds in the binding 
sites because at least one carboxylate was involved only in H 
bonding. This suggests that ionic bonds, despite their impor- 
tance, are not sufficient to stabilize citrate-camosinase com- 
plexes and that H bonds, van der Waals interactions, and mutual 
adaptability also play a key role. This model may also explain 
the high specificity of the citrate recognition sites in camosinase. 

Fluctuations in the Free versus Complexed Camosinase. 
Four MD simulations were performed involving camosinase 
alone, camosinase with three citrates, camosinase with car- 
nosine, and camosinase with citrates plus camosine. Fkst, die 
root-mean-square deviation (RMSD) profiles of camosinase 
backbone atoms was analyzed relative to the initial stmcture 
(i.e., after the heating phase). This was undertaken to assess 
fluctuations during simulations. Figure 4 shows that the four 
MD shnulations yielded similar RMSD profiles, tlius confimaing 
the overall stability of all systems. Although 5 ns may not be 
sufficient to entirely explore the conformational space of the 
structure of unliganded serum camosinase, its conformational 
rearrangements described below clearly appear to be mainly due 
to the dynamic response of the enzyme and not to random 
stmctural distortions. In particular, the sole RMSD profile, which 
is somewhat different, is that seen with camosinase in the sole 
presence of three citrates, which shows slightly higher RMSD 
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Figure 5, Percent change of each type of secondary sunicture in the 
three trajectories, compared to the trajectory of camosinase without 
the ligand (zero values). White bars: camosinase with citrates. Grey 
bars: camosinase with camosine. Black bars: camosinase with citrates 
plus camosine, 

values postulated to be due to m unfolding effect caused by 
the citrate ions. 

This effect of citrate ions was confirmed by separately 
calculating the percentage of each type of secondary motif. The 
last 4 ns of each trajectory was monitored by sampling 40 frames 
(1 frame/100 ps) and calculating the percentage of each type of 
secondary motif for each frame. The changes were expressed 
as differences between the free (taken as the zero value) and 
complexed camosinase (namely, camosinase with citrates, 
camosinase with camosine, and camosinase with camosine plus 
citrates). 

As shown in Figure 5, the simultaneous binding of camosine 
plus three citrates (black columns) caused an effect that was 
not the mere sum of the separate effects of citrates (empty 
colunms) and camosine (grey columns). Indeed, the citrates 
bound in the absence of camosine induced a partial unfolding 
of the protein with an increase in coll stractures, whereas the 
camosinase— camosine complex showed a modest decrease in 
the a-helix and turn motifs, a decrease compensated by an 
increase m strand and coil stmctures. 

In contrast, camosine and citrates together (black columns 
in Figure 5) induced an- increase in the percentage of a-helix 
and strand motifs, with a corresponding decrease in the 
percentage of bridge and turn motifs. In particular, the increase 
in a^helix abundance distinguishes this complex from the others, 
suggesting that the simultaneous binding of camosine and 
citrates induces a net shift from turn to a-helix motifs* 

To reveal which protein segments were affected most by the 
presence of citrates, camosine, or both, the relative time spent 
by each residue in a given secondary motif was also calculated 
over the 40 sampled frames. This was done for free camosinase 
and its three complexes. The results are expressed as the percent 
change relative to free camosinase (taken as the zero value). 
Figure 6A and B show that the simultaneous presence of 
camosine plus citrates induced a clear tendency from turn to 
a-helix rearrangements in the His390-Ala404 segment (Figure 
6A and B), which is a connecdng region between the lid domain 
and the catalytic domain (Figure 6C). 

Moreover, Figure 7A .shows that carnosine and citrate to- 
gether induced a marked stabilization of the a-helix stmcture 
of the His354-Ser373 segment,- which is located in the upper 
surface of the lid domam (Figure 7B). The stabilization of the 
lid doinain and the changes in the region connectmg the two 
domains are postulated to be implicated in a confomaational 
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Figure 6. Conformational changes observed during the MD simulations in the His390-AIa404 segment during the diree trajectories. The results 
are expressed as the percent of change relative to free camosinase (taken as the zero value). Whire bars: camosinase with citrates. Grey bars: 
camosinase v;^ith camosine. Black bars: camosinase with citrates plus camosine. (A) Percent changes in relative time spent in an a-helix envixonmeat 
by each residue in the His390— Ala404 segment* (B) Percent changes in relative time spent in a )8-tum environment by each residue in the His390- 
Ala404 segment. (C) Position of the Hts390— Ala404 segment in die connecting region between the two domains of camosinase. 



shift, which can modulate the accessibility of the catalytic site 
(see belov^). 

The changes in secondary structure around all carnosinase 
residues were carefully checked, but the results were not as clear 
as those described in Figures 6 and 7, and their analysis did 
not reveal any other significant conformational transition. In 
other words, the conformational changes in the considered 
segments demonstrate that a small molecule such as die citrate 
ion can induce significant and specific stractural e&cts in 
protein regions far away from its binding sites. 

Evidence for an Allosteric Effect of Citrate Ions. Interest- 
ingly, the most evident conformational shifts observed in the 
camosinase—camosine-'citrate complex concern the lid domain, 
whereas no marked conformational changes were seen in the 
catalytic domain. To shed light on the possible effects of citrate 
on the catalytic domain, die docking scores for the camosinase-' 
camosine complexes with and witliout citrates were evaluated 
by sampling one frame each for 0.1 ns for the first ns, one jOrame 
each for 0.4 ns for the remaining 4 ns, and submitting them to 
FlexX to obtain the corresponding scores. 



Figure 8 shows the score profiles as calculated during die 
two simulations, revealing that the presence of the three citrates 
in their respective binding sites leads to better scores for the 
binding of camosine. Indeed, in both cases, a worsening of 
docking scores was observed during the first ns (the equilibration 
phase) mainly due to an initial rearrangement of the catalytic 
site. But in the presence of citrates, this worsening was overcome 
by a progressive improvement of score values that reached their* 
minimum after about 3 ns and dien remained lower by about 
20 kcal/mol than the scores without citrates. Noticeably, in both 
simulations, the best score was reached after 3 ns, but the 
difference between the two lowest scores (—29.99 kcal/mol 
without ciu-ates vs "54.77 kcal/mol in the presence of citrates) 
affords impressive evidence for an effect of citrates on die- 
affinity of camosinase for camosine: 

When analyzing the various complexes and dieir stabiUzing 
interactions, we found that the presence of citrate ions induced 
a Significant increase in hydrophobic/van der Waals contacts. 
This may mean that citrate ions are able to influence the catalytic 
site by unproving the mumal adaptability between the enzyme 
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Figure 7. Conformational changes observed during the MD simulations 
in tlie Hi5354-Ser373 segment during the tliree trajectories. The results 
are expressed as the percent of change relative to free camosinase (taken 
as the zero value). White bars: cai'nosinase wjtb citrates. Grey bars: 
camosinase with carnosine. Black bars: carnosina^e with citrates plus 
cmosine. (A) Percent changes in relative time spent in an a-helix 
environment by each residue in the His354— Ala373 segment. (B) 
Position of die His354-Ser373 segment in the upper surface of the Hd 
domain, 

and substrates. Clearly, the key polar interactions that activate 
camosinase (e.g., the interaction between zinc ions and tlie 
carbonyl group in the substrate) are independent of the presence 
of citrate, but the increased carnosinase affinity induced by 
citrates can be explained by the emergence of several additional 
hydrophobicA^an de Waals interactions. 

To validate the results with citrates, we performed 1 ns MD 
simulations of the camosinase-- carnosine complex in the 
presence of (1) only one citrate molecule (one simulation for 
each citrate position) and (2) 9 chlorine ions in place of the 3 
X 3 carboxylate groups in the citrates. The results indicated 
very modest effects on the folding of serum camosinase (data 
not shown), suggesting tliat (1) the effects of citrates cannot be 
explained by mere electrostatic shielding effects and that (2) 
all three citrates are simultaneously required to obtain the 
described allosteric effect This last finding is of interest because 




Figure 8. Effects of the presence of three citrates on the docking score 
of die carnosinase— carnosine complex during the molecular dynamics 
simulation. The graph reports the docking scores calculated firom tlic 
frames stored everyO.l ns during the first ns and every 0.4 ns during 
the remaining 4 ns- 

it suggests a synergetic cooperation between three citrate ligands 
and invites experimental verification. 

Comcluslotis and Outlook 

In our opinion, the most unexpected information to emerge 
from this study is the molecular mechanism by which citrate 
ions act as allosteric activators of human serum camosinase. 
Allosteric transitions were discovered in the early 1960s by 
Monod and his eminent colleagues.^»-^ Since then, innumerable 
articles have been published showing the phenomenon to exist 
in enzymes, receptors, transporters^ and ion sensors.^'^^^ Tools 
applied to assess allosteric transitions include X-ray crystal- 
lography, enzyme kinetics and functional responses, protem 
engineering, and site-directed mutagenesis. To the best of our 
knowledge, little if anything has been done in silico using 
molecular modeling and molecular dynamics approaches. This 
article could therefore be one of the first reports detailing the 
molecular mechanism of an allosteric enzyme activator using 
molecular dynamics simulations. 

Inhibition of serum camosinase may be a useful therapeutic 
approach in the treatment of diabetic nephropathy. Recently, it 
was found that a trinucleotide repeat in exon 2 of the serum 
camosinase gene CNDPl ,was associated with nephropathy, 
whereas the shortest allelic form (CNDPJ Mannheim) associated 
with lower serum camosinase levels was more conomon in the 
absence of nephropathy.^^ In other words, the number of leucine 
repeats in the leader peptide of the CNDPl semm camosinase 
gene appears to be associated with susceptibility for diabetic 
nephropathy. These results also suggest that tlie design of semm 
camosinase inhibitors could afford renal protection in diabetes, 
especially, in the most exposed patients having an increased 
expression of camosinase due to the trinucleotide repeat in the 
CNDPl exon 2. Our homology model of human semm 
camosinase and especially of its catalytic site appears specific 
and explicit enough to allow the design of competitive inhibitors. 

Furdiermore, carnosine has been found as a bioactive peptide 
in several patliological conditions, but its rapid hydrolysis in 
blood limits its pharmacological application as a detoxifying 
agent toward cytotoxic metabolites involved in tlie pathogenic 
mechanisms of atherosclerosis, ischemia, and nepliropathy. It 
would therefore be of potential interest. to design camosine 
mimetics that retain the favorable effects of- camosine but are 
resistant to camosinase. Again, our homology model offers a 
first tool toward this objective. 

There is tlius a large body of evidence indicating that serum 
camosinase can be considered a promising and novel dmg target 
' and that therapeutic effects can be achieved either by inhibiting 
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or activating the enzyme depending on the pathological state. 
Activation of serum camosinase might prove an approach to 
treat camosinemia/homocaniosineinia, a very rare (prevalence 
1:500.000) inherited and incurable metabolic disorder character- 
ized by severely impaired neurological functions and develop- 
mental delays, caused by seram camosinase defidency and 
leading to high serum and urine levels of homocamosinc/ 
camosine. Dietary restriction of carnosine was found to reduce 
homocamosine and camosine to normal levels but was inef- 
fective or only slightly effective in improving the general state 
of patients, suggesting that neurological disorders are not 
associated with high serum levels of hisiidine-containing 
dipeptides but may be caused by an impaired metabolism of 
dietary camosine thai probably affects the content of histidine- 
containing dipeptides in tissues. Hence, a promising therapeutic 
approach might be to partially restore the activity of serum 
camosinase by an as yet undefined mechanism. The results 
presented here bring consistent evidence for an allosteric 
activation of serum camosinase, identify three allosteric sites 
in the protein, and suggest a molecular mechanism by which 
the enzyme's activity is increased. 

But whether inhibitors or allosteric activators, the design of 
coimpounds acting on camosinase would be of significant value 
to gain a better understanding of the physiopathological roles 
of the enzyme and the substrates whose concentration it 
regulates. Indeed, the biological function of these substrates 
remains poorly known, and findings of clinical importance may 
be down the road. For example, a significant decrease in serum 
camosinase activity has been found in patients undergoing heart 
surgery, suggesting tiiat this might be a functional mechanism 
to protect against ischemic brain damage.^^ 

Computational Methods 

Construction of Camosinase. The amino acid sequence of 
human serum camosinase (CNl) was retrieved from the Genebank/ 
EBI Data Bank (accession number: AX1397473^). The primary 
transcription sequence was submitted to the Signal IP 3.0 Server^^-^^ 
to predict die cleavage site of the signal peptide. Predictions show 
the highest probability of cleavage between Ser27 and Ser28. 

The primary sequence was submitted to Fugue,^^ an on-line 
sequence-structure homology recognition software that generates 
realistic models recognizing distant homologies by sequence- 
structure comparison. Tlie best model produced by Fugue was based 
on the structure of PepV dipeptidase from Lactobacillus del- 
hmeckiP (pdb id: ILFW) as the template. Because this protein 
has been cocrystallized with an inhibitor (namely, AEP, 3-[(l~ 
amino-a-carboxy-ethyO-hydroxy-phosphinoyll-i-methyl-propion- 
ic acid), the catalytic site was poorly accessible and could not 
accommodate any ligand without distorting its confomaation. Hence, 
this model was discarded, and a more suitable model was searched 
for among the 20 solutions proposed by Fugue. 

Results from the second to seventh rank were ignored because 
they predict less than 40% of the camosinase sequence, and thus, 
the choice was for the eighth model based on the experirnenial 
structure of /?-alanine synthetase from Saccaromyces kiuyverP (pdb 
id: IR3N) as the template. This choice was justified by three main 
reasons: (1) tiiis model predicts die highest percentage of camo- 
sinase sequence leaving a few shoit unpredicted fragments; (2) the 
homology percentage between human semm camosinase and 
/3-alanine synthetase (as computed by using pairwise alignment) is 
significant (52.6%) and comparable to that with PepV dipeptidase 
(58.3%); and (3) camosine contains )3-alanine in its stmcture, and 
thus, some structural analogies between the catalytic site of die 
enzymes (namely, camosinase and ^-alanine. synthetase) are prob- 
able. What is more, tlie model based on j8-al*anine syndietase shows 
a more accessible catalytic site, and indeed, docking analyses and 
MD simulations described here are performed with this model. 
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Specifically, the missing fragments in die selected caroosinase 
model were Serl-Ala?, Ala67-Asp78, ArglU-GlyllS, Argl90- 
lle200, Val37J-Asn379, and Phe431-lie434. The four longest 
fragments (i.e., Serl-AlaT, Ala67-Asp7$» Argl9Q-Ile200, and 
Val371-Asn379) were separately resubmitted to Fugue, selecting 
the model best fitting the remaining backbone, whereas die two 
shortest fragments (i.e., Argll4-Glyl 15 and Phe431-Ile434) were 
manually added using the VEGA software.^^ 

Widi the full-length backbone constructed, side chains and 
hydrogen atoms were added using VEGA. To remain compatible 
with physiological pH values, die side chains of Arg, Lys, GIti, 
and Asp were ionized, whereas His residues were considered neutral 
by default The complete model was carefully checked to avoid 
unphysical occurrences, such as cis peptide bonds, wrong configu- 
rations, or colliding side chains. In these first phases and in the 
following steps, die model soundness was assessed using Procheck*® 
and VerifySD."*^ A preliminary energy minimization was performed 
on the complete model to avoid high-energy interactions. During 
this minimization, the backbone atoms with the exception of atoms 
in the added fragments were kept fixed to preserve global folding. 

To correctly arrange the zinc ions in die catalytic site of die 
model, metal ions and ligating residues were extracted from the 
stmcture of /^-alanine synthetase and inserted in the camosinase 
model by superimposing die Ca carbon atoms of die coordinating 
residues to die corresponding atoms of camosinase. All of these 
atoms except for the zinc ions were dien deleted. Finally, an energy 
minimization was performed by fixing the atoms outside a sphere 
of 12 A radius, centered on die middle point of die segment 
connecting die metal ions. 

To gain a better relaxation and a more correct arrangement of 
the complete serum camosinase model, a molecular dynamic 
simulation was performed in vacuo. The simulation was divided 
into a heating phase from 0 to 300 K (3 ps, i.e., 1 K/IO iterations) 
and an equilibration phase (0.5 ns) with constant temperature (300 
db 25 K). During the simulation, only the zinc ions were kept fixed, 
Tlie last frame of this simulation was further minimized and 
represented the starting point for the calculations to foDow. This 
optimization procedure yielded a refined structure, which differs 
from the starting one for a RMSD equal to L67 A (as computed 
from only the backbone atoms), suggesting an overall stability of 
the predicted folding of camosinase. 

Building of Ligands. A set of five histidine-containing dipeptides 
was considered (Chart 1). namely, carnosine (^-alanyl-L-histidine), 
homocamosine (y-aminobutyryl-L-histidine), carcinine (/J-alanyl- 
histamine), A^-acetylcarnosine,*and anserine 05-alanyl-3-methyl-L- 
histidine). All molecules were built using the VEGA software.^^ 
Carboxylic and amino groups were ioni2ed, whereas the histidine 
imidazole ring was neutral (N*'-H tautomer). After a preliminary 
mmhnization, partial charge attribution and geometry optimization 
were perFonned by die MOP AC 6.0 program (keywords = AMI, 
PRECISE, GEO-OK, MMOK).^^ Tlic ligands underwent a Monte 
Carlo confomiational seai-ch using the Quanta/CHARMm package 
(MSI, Burlington, MA), which generated 1000 conformers by 
randomly rotating the rotors. All conformers were minimized to 
avoid high-energy interactions, and the lowest energy one was used 
in docking calculations. The citrate ion was built in its tricarboxylate 
form using the same procedure. 

Docking of Ligands. The software PlexX was used lo perform 
docking analyses."" FlexX is a fast-automated docking program that 
takes ligand conformational flexibility into account by an incre- 
mental fragment technique. Docking analyses involved residues 
enclosed in a 12 A radius sphere centered on the metal ions, 
generating 30 complexes for each ligand. Even if the exact 
mechamsm of hydrolysis can only be postulatcd,^^ the docking 
parameters can be justified considering that at least one metal ion 
must be involved in catalysis, and thus, it is 'highly probable that 
substrates- bind close to the zinc ions, 

.Docking .of Citrate. To investigate the role of citrate as an 
allosteric activator, the tricarboxylate form of citric acid was docked 
into the camosinase model. The docking was performed using 
FlexX^^ and involved the entire protein stmcmre, but the metal ions 
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were removed to avoid an overestimation of electrostatic interactions 
that would have concentrated the best docking solutions close to 
the ions. As with substrate docMng, 30 poses were generated. 

Molecular Dynamics. Molecular dynamics simulations were 
perfonned with the following models: (i) camosinase alone; (2) 
caraosinase with.three citrate ions; (3) carnosinase with camosine; 
and (4) camosinase with camosine and three citrate ions. All models 
were surrounded by a 10 A layer of TIP3S water molecules as 
generated by the software VEGA.3S> Cylindrical boundary conditions 
(radius 45 A; height 124 A) were applied to the simulation space. 
The hydrated complexes underwent a preliminary minimization of 
the relative position of the solvent molecules to .eliminate any high- 
energy interactions. 

Molecular dynamic simulations were organized into an initial 
period of heating from 0 to 300 K over 3000 iterations (3 ps. i.e., 
I K/10 iterations) and a 5 ns simulation phase with constant 
temperature. Only frames memorized during the simulation phase 
were analyzed. To differentiate the allosteric effect of citrate from 
a mere charge-shielding effect, a MD simulation of the camosi- 
nase-camosine complex was carried out with 9 chloride ions 
instead of the 3 x 3 citrate caihoxylates. The simulation lasted I 
ns with the same computational details. To better explain the effects 
of citrates, three MD simulations of the camosinase'-carnosine 
complex with a single citrate (one MD for each citrate binding site 
being occupied) were performed. The simulations lasted 1 ns with 
the same computational detaOs. 

Computational DetaOs. The software NAMD''^ was used to 
perform energy minimizations and molecular dynamics, using 
Gasteiger's atomic charges and the atom types in the CHARMM 
v22 force field. The metal ions were kept fixed in all simulations. 
The MD simulations were analyzed using STRIDE^^ to determine 
the relative proportion of secondary structure motifs in the various 
frames. 

All simulations had the following characteristics: minimizations 
with the conjugate gradients algorithm, convergence limit (RMS) 
= 0.01, maximal number of iterations = 5000; molecular dynamics 
with constant temperature in the range 300 ib 25 K by means of a 
thermostatic bath according to the Langevin's algorithm; and 
Integration of Newton's equation each fs according to Verlet's 
algorithm. The frame stored 5000 iterations (5.0 ps) each, yielding 
1000 frames per trajectory. 
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